MATHEMATICAL ANALYSIS SECTION 
MAS REPORT NO. 1 



'THE TRW- 130 (AN/UYK-i) APPLIED TO 
TRACKING RADAR DATA PROCESSING" 



Loren D. Enochson 



15 August 1962 



THOMPSON RAMO WOOLDRIDGE INC. 

RW Division 

CANOGA PARK, CALIFORNIA 



CONTENTS 



5. 



Introduction 

1. 1 Problem Statement 

General Description of the Atlantic Missile Range 
Acquisition System (ARCAS) 

Detailed ARCAS Program Description 

3. 1 General 

3. 2 Input-Output and Real Time Orientation 

3. 3 Designation Mode 

3. 4 Acquisition Mode 
Mathematical Details 

4. i Coordinate Conversions 

4. 2 Affine (Site to Site) Transformation 

4. 3 Position Updating and Data Validity Checks 

4. 4 Refraction Correction 

4. 5 Scan Mode Selection 

4. 6 Coordinate System Selection 
Program Extensions 

5. 1 Smoothing 

5. 2 Orbital Parameter Calculations 
References 



Page 

1 

1 



1 

3 

3 

4 

5 

6 

7 

7 

8 

12 

14 

15 

16 

17 

17 

20 

22 



A major portion of the material included in Section 3 was prepared by 
Isabelle Wengert, Information Systems Department, Thompson Ramo 
Wooldridge Inc. , R"W Division. 



THE TRW- 130 (AN/UYK-1) APPLIED TO 
TRACKING RADAR DATA PROCESSING 

1. INTRODUCTION 

A program presently exists to operate AN/UYK-1 computers which 
assist radars in the Atlantic Missile Range in acquisition and tracking of 
a target. This program is described along with mathematical details 
necessary for the implementation of such a program. Also, possible areas 
of sophistication are indicated to improve the effectiveness of such a program 
in some respects. These modifications, however, require increased com- 
puting time in general. 

1. 1 Problem Statement 

The basic problem is to automate those functions necessary for a 
tracking radar to acquire a target and then begin automatic tracking of 
this target. It is desirable for this target acquisition to occur as soon as 
possible in order that the radar will be able to output good position data 
which may then be used by a large central computer in such areas as the 
range safety functions (i.e. , impact prediction). The application of a 
digital computer to this problem allows a high speed (real time), accurate 
solution. 

Several features of the AN/UYK-i make it highly suited to this 
problem. First, it is a relatively high speed computer which, along with 
the adaptability inherent in the stored-logic feature, allows the necessary 
functions to be performed in the short cycle time (100 ms) allowed. 
Secondly, it is "ruggedized" which, due to the relative remoteness of the 
sites, eases the associated transportation and environmental problems. 
Thirdly, the AN/UYK-1 has a sophisticated interrupt system which allows 
response to the associated external devices in synchronism with real time 
with a minimum of programming problems. 

2. GENERAL DESCRIPTION OF THE ATLANTIC MISSILE RANGE 
ACQUISITION SYSTEM (ARCAS) 

The basic make-up of the system is: 

1. Six Radar-Buffer-AN/UYK-1 Computer sites on the Atlantic 
Missile-Range. 

2. A submarine cable for data transmission from' site-to- site. 

3. A high-speed data transmission system from radars to a 
larger computer at the Cape Canaveral site. This computer 

will transmit data over the submarine cable just as an AN/UYK-1. 

The general objective is to be able to have a site begin tracking a missile, 
obtain position data from this radar site, and then transmit this data down- 
range. This data must be presented to a downrange radar such that it will 
be properly directed at the missile in order to allow it to lock-on the target 
and go into an automatic tracking mode. 



Also, the position data obtained from the various site radars will be 
sent directly to the central computer which will use the data for impact 
prediction and associated range safety functions in addition to transmitting 
position data as an AN/UYK-1 does. Also, the central computer and each 
AN/UYK-i must present data suitable for plotting. Namely, in cartesian 
(X, Y, Z) coordinates referenced to the local radar site. 

The AN/UYK-1 computer fulfills these objectives by performing the 
following functions: 

1. Read range, azimuth, and elevation data (RAE) from a buffer 
which is connected between the local radar and the computer. 

2. Test this data for reasonableness. 

3. Apply certain corrections to the data. 

4. Convert it to a cartesian form' referenced to the local site. 

5. Update this data for various time delays in the system. 

6. Determine which one of five coordinate system origins is 
approximately closest to the target. 

7. Transform the data to this common reference origin. 

8. Insert a set of ID bits to specify this origin. 

9. Scale the data to appropriate units and round to 12 bits. 

10. Output the data along with a scale code and origin code to a 
buffer which will transmit it to the submarine cable. 

Items (1) through (10) are called the designation or tracking mode. 

11. Read cartesian (X, Y, Z) data from a buffer which is referenced 
to some common origin. 

12. Scale the data based on scaling code bits. 

13. Determine this origin from the ID bits. 

14. Select the proper matrix to transform this data to the local site. 

15. Perform a reasonableness check on the data. 

16. Output this local X, Y, Z data to a buffer for plotting purposes. 

17. Convert this cartesian data to polar form (RAE). 

18. Update this data for time delays. 

19. Determine a scan mode recommendation. 

20. Apply corrections to the data. 

21.. Output this R, A, E data along with scan recommendation bits to 
a buffer which will give the data to the radar. 

Items (11) through (21) are known as, the acquisition mode. 
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The requirements that must be satisfied while performing these 
functions are: 

1. The entire program, including both modes, must be able to be 
executed in less than 100 milliseconds. 

2. The X, Y, Z data from the designation mode must be output to 
the buffer prior to the 50 millisecond point. 

3. The entire program must be contained in 5460 cells or less. 

4. A provision for an alternate mode called the fixed mode must 
exist which would result in transforming data to an origin 
referenced at the local site rather than a common range origin. 

3. DETAILED AR CAS PROGRAM DESCRIPTION 

3. 1 General 

Presented in this section is the basic program description with many- 
mathematical details omitted. Reference may be made to Section 4 as 
needed by the reader to clarify points not precisely described in the follow- 
ing discussion of the program. 

Every 100 milliseconds, cartesian coordinates referenced to the 
common coordinate system with a code indicating which station is in the 
best position for tracking the object at a given time are output to the other 
five units. From this data with the associated identification (ID) code, 
one station can begin tracking, and output to the radar, polar coordinates 
in the local coordinate system with a scan mode recommendation. When 
radar data is read, it contains a code bit, designated CS, which determines 
whether the local transformation matrix is to be used. This bit is controlled 
via a switch on the radar operator's console. 

Each station has an assigned ID and matrix elements which are desig- 
nated as the local coordinate transformation for that installation. For 
example, Patrick Air Force Base might have an ID of 4. Then matrix 4 
and vector 4 would be stored as the local matrix transformation at the 
Patrick installation, while Cape Canaveral might have an ID of 3; hence, 
matrix 3, vector 3 would be the local matrix transformation at the 
Canaveral Station. 

An OT bit code, "on track, "is input with channel B radar data. 
When the radar is on track, no scan mode recommendation is made. Other- 
wise, the program recommends a spiral or one of two raster scan modes. 
The rate of change of azimuth and elevation angles determine whether the 
spiral or raster scan mode is recommended. The spiral mode traverses 
equal azimuth and elevation angles and begins searching in the center and 
spirals outward. The raster scan mode, a slower process, begins search- 
ing at one corner of the pattern. The two raster scan modes differ only in 
the rate of change of azimuth. The method for scanning is only allowed 
to be changed every minute. 



There are two starting points for the operational program. One 
sets up one set of matrix elements to rotate the X and Y coordinates to 
an earth centered vertical orientation for plotter output. The second 
entry sets up an identity transformation to output X and Y oriented to the 
radar tangent plane. 

The entire program is executed in a maximum of about 85 ms and 
requires about 4200 cells of storage. Operations are performed in either 
single or double precision as dictated by the precision of the input or output 
data. The input data is, in the case of the AN/FPS-16 radar, 20 bits, 19 bits, 
and, 19 bits for range, azimuth and elevation respectively. This requires 
double precision operations to be performed in most of the designation 
mode operations such as the polar to cartesian coordinate conversion. The 
final X, Y, Z output to be transmitted via the submarine cable, however, is 

12 bits and a sign bit which allows certain scaling and output operations to 
be performed in single precision. The input to the acquisition mode is 
therefore 12 bits and a signfor X, Y, Z while the output for the radar is 15 bits, 

13 bits, and 1 3 Jbit's' for -•RrA, and,E respectively. This allows many 
functions of the acquisition mode, such as computation of A and E from 
X, Y, Z, to be performed in single precision. Some double precision 
arithmetic is required, however, and at times the necessary scaling of 
the numbers makes double precision operations convenient although not 
absolutely necessary. 

3. 2 Input-Output and Real Time Orientation 

Channel A inputs, through the buffers, the X, Y, Z data and the ID 
and SC codes, all needed in the acquisition mode. Channel A output 
consists of the computed X, Y, Z (both for transmission and for plotting) 
and the ID and scale code (SC) determined in the designation mode. 

Channel B inputs the R, A, E data and the CS and OT codes used in 
the designation mode. The output on channel B is R, A, E and the scan 
mode (SM) recommendation which is calculated in the acquisition mode. 

All input and output is performed in response to interrupts from the 
buffering devices and full use is made of the interrupt capabilities of the 
AN/UYK-i. This interrupt system also allows a limited self checking 
operation to be performed on the system. That is, a check may be 
performed at the conclusion of a program cycle to determine if the correct 
number of interrupts were received from each channel. If not, the 
results of that cycle may be incorrect, and synchronism with real time 
may have been lost. However, instead of a collapse of the complete 
system, if this was merely an intermittent failure, the AN/UYK-1 may 
reorient itself in the real time cycle and correct operation is recovered 
with possible loss of only a few cycles of data. 

This real time orientation is also the first step performed by the 
program when the system is first started. The orientation is accomplished 
by checking for the existence of input and output data requests (which 
remain in effect once initiated by the buffers), counting elapsed time to 
the next requests, and finally determining a channel A input data request 
which begins the program cycle. 



3. 3 Designation Mode 

Data Input 

Channel A data is input and the codes are extracted. Data is scaled 
by a power of 10, according to the scale code (SC) read with the data, and 
the units are changed for internal usage. Channel B data is input and the 
codes are removed. Data is stored, appropriately scaled, , as specified by 
the program. 

Validity Test and Extrapolation 

A counter, which is incremented every cycle, indicates whether 
enough points have been input for validity checking and for extrapolation. 
These functions are bypassed until the eleventh cycle. 

An average first difference is computed over the five most recent 
data points. This increment is then used for a single cycle extrapolation 
to use as a check value in the validity test, and also for time delay extrapo- 
lations. 

If a piece of data fails the validity test, an error indicator is set, 
and the check value is employed in subsequent computations for that cycle. 
However, the next two data points are not checked and are accepted for 
use in computations. The third point is again checked. 

Refraction and Bias Correction 

The program next adds a refraction correction, computed during 
the previous cycle in the designation mode. Constant bias corrections 
are added to R, A, and E. 

Polar to Cartesian Conversion 



The sine and cosine of azimuth and elevation are computed and then 
used to compute X, Y, Z. The intermediate product ALPHA = R cos E 
is stored for the refraction correction calculations in the designation mode. 
X, Y, Z are stored with the same binary scaling as range. 

Coordinate Transformation to Common 

If the coordinate system bit (CS), which is read in with channel B 
data, is zero the local coordinate transformation is used. Otherwise, 
the reference origin counter is tested. If the counter is negative, then 
X-Y is compared against four barriers. The barriers represent imaginary 
planes dividing the tracking range of one station from the next one. When 
X-Y is less than bar I, Matrix I is used in the transformation, I is stored 
in identification bits (ID), and a program transfer is set so that the program 
automatically uses the matrix I transformation for the next 9 cycles. After 
completion of 10 cycles, X-Y is again compared against the barriers and 
a new I is found. 



Maximum Coordinates 

Following the coordinate system transformation, the program 
determines whether X, Y, or Z is the largest coordinate. The largest 
coordinate is compared against stored constants for finding a scale factor 
to reduce X, Y, Z to one word format. The scale code is stored for output. 

Channel A Output 

The scaled data has the ID and scale code added and is formatted 
and output on channel A. 

3. 4 Acquisition Mode 

Coordinate Transformation to Local System 

According to the ID input with the channel A data, a set of matrix 
and vector elements is selected for the transformation from the common 
to the local coordinate system. The data also has an additional trans- 
formation applied for proper orientation of X and Y for the plotter. These 
two values are then formatted and stored for output to the plotter at a 
later time. 

Validity Check 

A validity check is performed in the same manner as for the R, A, E 
data in the designation mode. A slight difference occurs since the expected 
range of random errors in the R, A, E data is known, but the random errors 
in the X, Y, Z are a function of the geometry of the situation and the allowed 
error for X, Y, Z is therefore based on the computed first difference. The 
procedure is explained in greater detail in Section 4. 3. The time delay 
extrapolations must be delayed until R, A, and E are available since different 
delays are compensated for in R than in A and E. 

Cartesian to Polar Conversion 

Range is computed from the square root of the squares of X, Y, and 
Z. Azimuth is the arctangent of the quotient of X over Y. Elevation is 
the arcsine of Z divided by R. The quadrant for azimuth is determined 
and the two angles are converted from radians to units of fractions of a 
circle used for channel B input. 

Refraction Correction Determination and Bias Corrections 

In the refraction correction calculations, cotE is computed from 
R cos E/Z, where R cos E is available from the designation mode calcu- 
lations. Next, a table look-up through thirteen values is performed with 
Z as an argument to find a factor (N - N)(A) where A is a scale constant. 

See Section 4. 4 for details. Finally, AE = (N - N) cotE which is then 

s 

added to E. Constant bias corrections are next added to R, A, and E. 



Extrapolation 

R, A, and E first differences are now computed and used for extrapo- 
lation. A and E are extrapolated an additional cycle since the acquisition 
servo units for A and E work as a function of their error from the desired 
position. 

Scan Mode Recommendations 

Provided that on-track bit is not set and that time delay counter 
shows that one minute or more have elapsed since the last scan mode 
recommendation, a new scan mode is selected. The average differences 
for azimuth compared against the average differences for elevation and 
against a predetermined constant provides the method for making a raster 
scan mode or spiral scan mode recommendation. That is, the relative 
rates of change of elevation and azimuth are the selection criteria. 

Data Output 

Range units are converted for output. The two angles are changed 
to one word format and scan mode hits are added to the data before 
formatting for output to the buffer. The local plotter data is output 
followed by the acquisition data. 



4. MATHEMATICAL DETAILS 

4. 1 Coordinate Conversions 

Radar to Cartesian Coordinates (RAE — *"XYZ) 



Assume azimuth (A) is measured from north, positive in a clock- 
wise direction. It is assumed that if in some unusual situation a target 
passes directly overhead such that the elevation angle (E) exceeds 90 
degrees, then the azimuth becomes 180 degrees out of phase. That is, 
the radar remains stationary in azimuth while the target passes directly 
overhead. In other words, there is a discontinuity (a singular point) at 
E = 90 degrees. The equations to convert from radar to cartesian 
coordinates then are: 

x = R cos E sin A 

y = R cos E cos A (4. 1) 

z = R sin E 

The change in the sign of cos E when E exceeds 90 degrees accounts for 
the effect of A becoming 180 degrees out of phase. 

Cartesian to Radar Coordinates (XYZ — HEtAE) 
The equations for this conversion are: 



2 2 2 

R = yx+y+z 

E = arcsin (z/R) (4. 2) 

A 1 = arctan (y/x) 
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The proper quadrant for A is determined by a check of the signs of x and y: 

x,y 



+> + 


A = 90° - A' 


-,+ 


A = 90°+ A' 


— 9 — 


A = 270° - A' 


+, - 


A = 270° + A' 



4. 2 Affine (Site to Site) Transformation 

Notation: Rotations about the various axes are as follows: 
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R = I cos sin 
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R Bz = 1 sin e cos e ° 



The affine transformation from site i to site j is denoted by 



(x). = A.'.(x). + T.. (4. 3) 

v 'j ij i ij 



where 




and A., is the rotation matrix with T.. being the translation vector. A 

y y 

position vector with a prime indicates earth centered vertical orientation 
and is related to (x) by 



Also define: 



(x): =R , (x) (4.4) 

i " °i* 1 

<p. = latitude of site i 
^ . = longitude of site i 

Alt- = X _ \ 
i J 

s. = distance to the equatorial plane along 
the normal to the tangent plane at 
site i. 
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d. = distance to center of earth along equatorial plane 
from intersection of tangent plane normal with 
equatorial plane 

h. = height above sea level of site i 

6. = angle between tangent plane normal and earth 
centered vertical axis at site i 



See Figure 1 for an illustration of these quantities. 




Figure i. Site to Site Coordinate Transformation 



Affine Transformation: The breakdown of the transformation from a 
tangent plane oriented system to an earth centered vertical system is 



(x). 1 = R r 
J -°i 



x 



R. 



6.x 
Y 3 



Afcy 



■<i.x 



{x) i +- 







0i..+ s/ J 
i i 




(4.5) 



This is the transformation used in the designation mode in going from 
the radar oriented cartesian system to the earth centered vertical 
orientation in which data is transmitted via the submarine cable. 



Let 



cAX 





-sAX 


- s0. sAX 


C0. 


-s0. cAX. 


co. sAX. 


s0. 


c^.cAX, 



(4.6) 



To get Eq. (4. 5) in the form of Eq. (4. 3), one combines terras and obtains 

(° 

(x).' = R . R^ R A \ R ^ (x). + R . R^ R A , R . 

x 'j -o.x 0.x AA-y -0.x i -o.x 0.x AA-y -0.x I 

J J x J J * \h. + s. 

v i i' 

n n / ° 

+ R -5.x R 0.x R AXy ° - R -5.x%.x ° " R -5.x ° 

'J J \d./ J ^ \d./ J Vh. + s./ 

Then 

A.. = R^ R A . R , (4. 7) 

ij 0.x A\y -0.x 

Letting c0 = cos and s0 = sin 0, the rotation matrix is 
cAX sAXs0. -sA\c0. 

A.. = I -s0. sAX C0.C0. + s0.cAXs0. c0.s0. - s0.cAXc0. j (4.8) 

C0. SAX S0. C0. - C0.CAXS0. S0.S0. + C0.CA\C0. 



B.. =R^_ x R AXy = j~s0 4 sAX c0, ~ S ; cAX| (4,9) 



The translation vector then becomes 

T.j = A ( j + B.. I J -R ' 0-1 | (4. 10) 

\h.+ s./ \d. / j \ d. / \h.+ s. 

ii i J J J 

The transformation from tangent plane orientation at site i to earth 
centered vertical orientation at site j used in the designation mode can 
now be written 

(x)l = R , A.,(x). + R , T„. (4. 11) 

x 'j -6.x ij N 'i -ox ij 

J J 
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All the matrices involved are rotations and therefore orthogonal which 
implies their inverses are given by their transpose. The inverse trans- 
formation from earth centered vertical orientation at site j to tangent 
plane orientation at site i for use in the acquisition mode is therefore 
given by 



(x), = 



R , A.. 
-o.x 11 

J ■ 



(x)l -Re T.. 
J - 5 j x iJ 



where the T superscript indicates the transpose. Simplifying 
(x) ; = 



R. A.. 
-5.x 11 

J . 



T T 

(x)! - A.. 1 T.. 



(4. 12) 



(4. 13) 



One other relation is useful in calculating the necessary matrices for 
the ARCAS program which relates (x). and (x)! . Using Eq. (4. 11) 

(x.) = R s (x)! = R R (R , A..(x). +.-R K T..) 
x y o.x* 'i o-.x* -o.x if 'i -o.x ir 

J J J J J J J J 



= A..(x). + T.. 

y i ij 



(4. 14) 



Solving for (x). 



(x). ■= A.. T f(x). - -T..1 



Thus 



(x)! = R R (x). = R - |a.. T (x). - A.. T T..l 
'i -5.x x 'i -o^L ij v 'j ij ij J 



(4. 15) 



Therefore, since 



and 



A.. = A.. 
Ji y 



T.. =--A.TT.. 

J 1 y y 



all necessary transformations can be calculated from A., for i=*i, 

ij J 



R c _ for all i, and T... for i=»j. 
i 



■6.x ij 
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Assumed Shape of the Earth : 

The site latitude, longitude, and height above sea level, <f>, \ , and 
h are given from survey data. The angle 6, and distances s and d are 
calculated from based on an assumed shape of the earth. The 
assumed parameters are 

a = semi major axis = 6 975 459. 5 yards 
b = semi minor axis= 6 951973. 1 yards 
a- b _ 1 
a " 297 

The deviation of the normal 6, is calculated from the relation 

e(i - e/2) sin 20. 0. 003361335 sin 20. 

tan 6. = _i = _i_ 

1 -(2e - e 2 ) sinV 1 - 0. 006722670 sin 2 ©. 
The distances s and d are calculated from 

s b 2 

s. = , „ „ T7T (4 „ 17) 



sin 0. (b + a cot 0.) 



2 2 

cot 0.(a - b ) 

d. = 



1 /u 2 j. Z , 2 ^^/ 2 
(b + a cot 0.) 



4 - 3 Position Updating and Data Validity Checks 

An average first difference over five points is computed for use 
in updating and validity checks. This is given by 



x 5 " x l 




(4. 16) 



(4. 18) 



Ax= — r- = 4/ (x i+i "^ < 4 - 19) 



This linear difference is simple to compute and an example of its 
expected accuracy is given in Ref. [ll , 'It is showricthere 

that any error contribution is compatible with the remainder of the 
system accuracy. 

Extrapolation is performed merely by multiplying Ax by the 
appropriate time delay constant t,. This is used in the case of both 
x, y, z and R, A, and E. 
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That is, for example 



x = x c + t , Ax 
e 5 d 



where x represents the extrapolated value of the parameter x. 



e 



Validity checking is performed in the same way except that the 
latest value, x-» is compared against a one-time -interval extrapolated 
value x , based on the previous five points x n » x,, ... , x . . For x, y, z 
the test is 

|x - xj* 3 Ax + € (4. 20) 

where e is some small numerical value employed to avoid comparisons 
against zero. The ± 3 Ax limits allow for extrapolation and random 
errors while rejecting errors due to lost bits and the like. 

Since basic expected random errors in R, A, and E are known, the 
standard deviation of these distributions may be used in the limits, and 
the validity check is 



|R v -R 5 |-4.5cr R (4.21) 

Similarly for A and E. 

The apparently wide 4. 5cr limits are chosen to minimize the 
probability of rejecting any legitimate data. The object of these validity 
tests is only to check for reasonable data. That is, one wants to reject 
data if a significant bit has been lost in the data transmission. However, 
no real data smoothing or anything of the like is intended in an operation 
Of this kind. 

An additional procedure is necessary in the validity check to avoid 
problems resulting from certain situations which may arise. If, somehow, 
a legitimate erratic change in course of the vehicle being tracked occurs, 
a point of data might be rejected. If this occurs, all data from that time 
on might be rejected and extrapolated values used. This could possibly 
take place when a missile is sitting on the launch pad and constant position 
data input is being received which causes the first differences to be zero. 
At lift off, the change in position could conceivably be significant enough 
to cause data to be rejected. To avoid this, if a piece of data is rejected, 
the subsequent two points are not checked but are accepted regardless. 
The validity check is then again applied to the third following point. 
Certain fluctuations might result in the computations over five or six 
cycles if two or three bad pieces of data were received in sequence. 
However, the probabilities of this are small if the error is just inter- 
mittent. The most probable situation of a single lost bit followed by good 
data is handled properly. If data is bad over a long period of time, the 
system has most likely failed and nothing can be done in the validity 
check to take care of this situation. 
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4. 4 Refraction Correction. 

The equation for use in the refraction correction is as follows and 
may be derived using Snell's law and assuming a stratified atmosphere 
along with some simplification obtained by neglecting factors of small 
effect. 

.-6 



AE = - [n s - N(z)l cot E(10 ) (4.22) 



The factor N is the surface modulus of the index of refraction and 
N(z) is the average of the modulus of the index of refraction from 
ground level to the height z of the target. For large altitudes, N(z) 
approaches zero. Some typical values for AE from Eq. (4. 22) are as 
follows: 

E = 12° R = 25 N. M. z * 31, 500 ft. 

Based on AMR data N ■** 365 is a typical figure and the average value of N 
up to a height of about z = 31,500 ft. is N«260. Using Eq. (4. 22), 
AE ** -. 5 mils which corresponds favorably to a correction of about 
-. 65 mils (6400 mils = 360 degrees) based on a typical atmosphere 
given in Ref t 2 . 

If E = 3°, R = 20 N. M. , then z « 6240 ft. and _N« 320. The 
correction then is AE «*-. 9 mils corresponding to AE«#-1.0 mils 
given in Ref. [2J. Other calculations yield similar results indicating 
maximum resulting errors on the order of ± . 2 mils and an rms error 
of somewhat less. 

The method for implementing this is based on the fact that the 
geometry of the tracking situation changes slowly as far as the refrac- 
tion correction is concerned. Therefore, corrections in the acquisition 
and designation modes may be the same, and using one computational 
result from the designation mode to calculate the correction in the 
acquisition mode will not significantly affect the validity of the refraction 
correction. 

The correction procedure is as follows. First calculate 

. „ R cos E 
cot E = 



where R cos E is available from the designation mode and z from the 
acquisition mode. Then, using z (height) as an argument, a table 
look-up will be performed on thirteen values to obtain a precomputed 
value 

k[N g - N(z)] 10 -6 
where k is an appropriate scaling constant and N is a function of z. 
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Then calculate 



AE = - (k [n s - N(z)] 10" 6 j(cot E) 



In the ARCAS program these calculations are most conveniently per- 
formed in the acquisition mode with the same AE used for the desig- 
nation mode. 

4. 5 Scan Mode Selection 

Typical scan modes available for a tracking radar are: 

I. A spiral scan out to 5 degrees from beamcenter. 

II. A raster scan covering 2.4 degrees in azimuth by 8 degrees 
in elevation with 0. 8 degree beam center spacing. 

III. A raster scan covering 0. 6 degree in azimuth by 8 degrees 
in elevation with 0. 6 degree beam center spacing. 

The decision to employ the radar spiral or raster scan mode can 
be based upon the rate of change of the target azimuth and elevation 
angles as seen from the acquiring radar. The radar spiral scan mode 
coverage differs from the radar raster scan mode coverage primarily 
by the fact that it subtends equal azimuth and elevation angles. An 

additional asset is that the radar beam begins searching at the center of 
the scanned area, where the target is most likely located, and spirals 
outward from the center. Therefore, the radar beam will, in all 
probability, acquire a target slightly sooner than with the raster scan 
mode which begins at one corner of the pattern. 

Proper selection of the radar scan mode can be automated by com = 
puting successive first differences of both the target azimuth and elevation 
acquisition look angles. The absolute magnitude of the elevation rate 
I El I can be divided by the absolute magnitude of the azimuth rate | Az | . 
If the result is less than 2. selection of radar scan mode I is automated. 
If the result is equal to or greater than 2. 0, then selection of radar scan 
mode II or mode III is automated as described below. The decision 
point of 2. is based upon the relative azimuth and elevation coverage 
ratio of 3. 3 for radar scan mode II. This decision point can be modified 
to reflect current operational practice if desired. 

The elevation angles scanned by the two raster modes are identical 
but the azimuth angles scanned differ by a ratio of 4 to 1. Therefore, 
the significant parameter is azimuth coverage. A raster scan mode 
decision based on this and the geometry of the situation can also be 
readily automated. Automatic computer selection of the best raster scan 
mode can be based upon the magnitude of the rate of change of azimuth 
angle | Az | . 

The rate of azimuth angle change that the AN/FPS-16 radar can 
accommodate, and still maintain reliable target track, is - 750 mils 
per second. Therefore, the radar raster scan mode can be 
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apportioned between 
coverage as follows. 



750 mils per second on the basis of azimuth 



Radar Azimuth 
Scan Coverage 
Mode of Scan Mode 



II 
III 



2.4' 
0.6 C 



Weighted 

Azimuth 

Coverage of 

Scan Mode 

0. 80 
0. 20 



Azimuth- 
Rate 
Coverage 

600 mils/second 
150 mils/second 



Based on these figures the following raster scan mode selection would 
be*. 

Range of Azimuth-Rates Scan Mode 
- 150 mils/second III 



150 - 750 mils/second 



II 



Note that the analysis employed here is not sophisticated and the 
problem of scan mode selection can be investigated in considerably 
greater depth. One can consider in much greater detail the actual 
beam coverage obtained with a certain scan mode, and the rate at which 
the scan is performed. Probability statements can then undoubtedly be 
made considering various possible situations which could arise in the 
geometry of the acquisition problem. The necessary analysis is by no 
means simple, however, and the selection criteria presented is satis- 
factory for most practical purposes. 



4. 6 



Coordinate System Selection 



The geographical layout of the AMR radar chain makes a coordinate 
system selection procedure easy to implement in real time. The selection 
of a coordinate system origin close to the vehicle being tracked is desir- 
able in order to minimize the magnitude of the numbers which are trans- 
mitted such that maximum precision is maintained. 

The coordinate origin selection is based on the fact that the radar 
chain lies on approximately a line of 135 degrees azimuth going down 
range. This approximation is sufficient for the origin selection purposes. 
Consequently, the barrier lines are 45 degree lines midway between the 
sites and have equations of the form 



x = y + b 



i = 1, 2, 3,4 



(4. 23) 



The b is the x intercept and the slope is unity. For a given site there 
are then four of these lines with four b intercept values. The procedure 
is then to compute |3 = x - y and compare with the 



b. # i = i, 2, 3, 4, 



Numbering down range, if 



b. 4 <(3<b. 
l-l — r — i" 



then the ith origin is selected. 
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5. PROGRAM EXTENSIONS 

The linear extrapolation technique described in Section 4. 3 is 
extremely well suited for real time calculations. However, future 
applications 'might demand greater accuracy for target acquisition and 
two areas of sophistication enter naturally. The first of these is data 
smoothing and the second is orbital type calculations. A short 
description of these two processes appears below. 

5. 1 Smoothing 

The process of data smoothing is performed in order to reduce the 
effect of random error (noise) in the given data. The mathematical 
procedure is that of fitting a polynomial (usually second degree) to the 
observations in the sense of least squares. The reason this procedure 
is suitable for real time applications is that a "smoothed" point may be 
represented as a linear combination of the observed points, with the 
coefficients being only a function of the number of points (if the time 
interval between points is constant). This same polynomial that is 
effectively fitted to the data may also be used for extrapolations and 
therefore the validity check may be constructed around the smoothing 
procedure. Smoothing is especially important in orbital parameter 
calculations where the most critical parameter is the velocity vector. 
Estimating the components of this vector requires as much error to be 
smoothed out as is practical. 

The derivation of the smoothing coefficients is presented below, 
Note that the analysis treats each coordinate separately which is not 
quite the same as a three dimensional fit accounting for all coordinates 
simultaneously. Also, the assumption is implicit that errors are from 
a stationary random process and uncorrelated. 

The assumption is made that a tracking observation x is of the 
form: 

2 

x i = a Q + a^t,. + a,t. + e . (4. 24) 

where e is a random error normally distributed with zero mean and 

2 
variance cr^ . That is, a second degree polynomial represents the 

data correctly, which is a valid assumption in practice. A smoothed 
point p. can then be given by: 

n 

(4. 25) 

where n + 1 is the number of points on which the smoothing is performed. 
A linear combination, which gives a best fit in the least squares sense of 
this type, is obtained as follows: Let n ..= 5 points for illustrative pur- 
poses. 
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The minimizing condition is 

4 

S = ) (x. - p.) = min (4. 26) 

fclJ 

The a., j = 0, 1, 2 are computed by taking partial derivatives and setting 
results equal to zero. That is, the so-called normal equations are: 

-ii_= 2S(x. - p.) = 

8 a i 1 



85 = 2St. (x. - p.) = (4.27) 



3aj 



iL. = 2St 2 (x. - p.) = 
8a 2 



Later calculations are simplified in the case of equally spaced time 
intervals when the points are relabeled. That is, let i = -2, -1, 0, 1, 2 
and upon "normalizing" time one gets t. = i, i = -2, -1,0, 1, 2. In 

matrix form the equations become 



n 


St 


st 2 


St 


st 2 


st 3 


st 2 


st 3 


st 4 




a o 



a 




a t | = | Ext ) (4. 28) 

2 



or more concisely 



Ta = x (4. 29) 



Note that by the judicious choice of indices the matrix T has zero ele- 

3 
ments. since St = and St = 0. Also, the coefficients a. will be 

J 
able to be expressed once and for all in terms of the observed x. since 

i f 
the t i are known. To be specific, if a., are the elements of T , 
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then 



2 2 



i=-2 i=-2 i^2 

2 



,. 2 

__ l 

i=-2 i=-2 i=-2 



l l = a 2i X Xi + a22 .Z ^ + Q23 Z ^ (4 " 30) 



a 2 = 



2- 2 



= a . \ x. + a_, \ x.t. + a.- VVx.t. 

31 L ' 32 ^ ii 33 ^ ii 

i=-2 i=-2 i=-2 



where 



_2 

i^2 

2 2 

i=-2 i=-2 



\ X. = (x 2 + x i + x n "•* X, + X_) 



^ x.t. = ^ x.i = (-2x_ 2 - x_j + x 4 + 2x 2 ) (4.31) 



A 2 v 2 - .2 



Zx.t. = \ x.i = (4x , + x , + x, + 4x _) 
li /_, i ^ -2 -i 1 -V 



i=-2 i=-2 



Assume a smoothed value at the midpoint, that is, for t- = 0, of 
•the given interval is desired. One sets t_ = in the polynomial 



P = a + a l t + a 2 fc < 4 ' 32 > 

from which one obtains 



p Q = (q u - 2a 12 + 4a 13 )x_ 2 + (a^ - a, 2 + a^Xj + Qj^ 
+ (a u + a 12 + a 13 )xj + (a u + 2a 12 + 4a 13 )x 2 



giving the smoothed value in the desired form. 



(4.33) 
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The numerical values of these coefficients are in this case 



b 1= 


-.08571430 


b 2 = 


. 34285714 


b 3 = 


.48571428 


b 4 = 


. 34285714 


b 5 = 


-. 08571430 



To extrapolate one interval, one only needs to set t = 3. The analysis 
for more points and for higher order equations proceeds in the same 
manner. A comprehensive discussion and an extensive tabulation of 
coefficients for many cases appears in Ref. 5 . 

Note that the implementation on a digital computer requires 
n + 1 multiplies and n adds for smoothing over n points. A highly 
improved procedure for real time use on a binary digital computer is 
described in Ref. [^ 3j and termed "almost least squares" smoothing. 
The procedure is to replace the b. by p\ given by 

p i . 

= *>i (4. 34) 

D 

where the p\ are a power of two and D is a normalizing factor. On a 
binary computer this replaces the multiply and add operations by shift 
and add resulting in a faster procedure. It has been shown (see Ref. [3]) 
that using n + 1 points in the almost least squares sense is about as 
efficient as n points in the least squares sense in terms of residual 
variance, 

5. 2 Orbital Parameter Calculations 

A further refinement in the acquisition problem is to implement 
position calculations based upon the parameters describing the motion 
of an orbiting vehicle. These are applicable, and highly desirable, in 
the case of tracking an orbiting satellite or when tracking a missile in 
the "free-fall" phase of its flights. Given the orbital parameters, the 
vehicle's position may be predicted with a high degree of accuracy for 
a given value of time by solving the equations for that value of time 
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even though the mathematical model is usually highly idealized. There- 
fore, for a sophisticated extrapolation and acquisition procedure, orbital 
calculations ■would be used. 

Assume a set of orbital parameters defined as follows, along with 
a given time t, is available: 

a ■* semi major axis of the ellipse 
e ■■» eccentricity of the ellipse 
cr - mean anomaly at epoch 
cos i - where i is the inclination angle 

!J2 - right ascension of the ascending node 
co - argument of the perigee 

Typical calculations to obtain position data in the form desired then are 
as follows. First, Kepler's equation must be solved for the eccentric 
anomaly E. That is 

E - e sin E = n(t + <r) (4. 35) 

where 



n =V5 



jjl = constant (4. 36) 



Note that an iterative procedure is necessary for obtaining E. The 
true anomaly f is then calculated from 



f = 2 tan" 1 



1 + e (E) 

tan - — - 



(4. 37) 



and the radial distance by 

r = a(l - e cos E) (4. 38) 

One now computes the coordinates in the orbital plane. 



x = r cos f 

y = r sin f 

The coordinate system for the orbital plane is then rotated to coincide 
with an inertial system. This is accomplished by a matrix multipli- 
cation. 



(4. 39) 
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That is: 



cos w cos ft - cos i sin ft sin w sin co cos ft + cos i sin ft cos 



CO 



x 



= cos co sin ft + cos i cos ft sin co cos i cos ft cos u - sin w sin ft ^ (4. 40) 

sin i sin co sin i cos co 

A further transformation of coordinates is necessary to reference 
this x, y, z to the local radar site followed by the conversion to radar 
coordinates. The above calculations along with those necessary for 
obtaining a set of orbital parameters from x, y, z position data may be 
found in Ref. [4]. 
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